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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TECHNICAL NOTE D-78O 

A METHOD FOR OBTAINING THE NONLINEAR AERODYNAMIC 
STABILITY CHARACTERISTICS OF BODIES OF 
REVOLUTION FROM FREE-FLIGHT TESTS 
By Donn B. Kirk 

SUMMARY 


A method is presented for obtaining the nonlinear aerodynamic 
stability characteristics of bodies of revolution from free-f light tests. 
The necessary conditions for the application of this method are (l) that 
the roll rate and damping encountered in a single cycle of oscillation 
be small , and (2) that the resulting motion be reasonably planar. Four 
approximations to the nonlinear restoring moment are considered and solu- 
tions are obtained in closed form: 

1. A single-term polynomial in an arbitrary power of the angle of 
attack 

2. A two-term polynomial having linear and cubic terms 

3. A three-term polynomial having linear, cubic, and quintic terms 

4. A three-term polynomial having linear, quadratic, and cubic terms 

An iteration procedure is formulated to allow the use of each of these 
approximations for obtaining the aerodynamic coefficients of bodies of 
revolution from free-flight test data. It is found that although the 
equations that are solved pertain strictly to planar motion, the solutions 
are applicable to motions that deviate to a fairly large degree from planar 
motion. 

Two of the approximations are applied to a set of data gathered from 
the Ames Supersonic Free-Flight Wind Tunnel. It is shown that one of these 
approximations is clearly superior to the other in fitting the basic data. 
The results of this better approximation are then compared with both data 
obtained from a wind tunnel and the results of two other methods for reduc- 
ing free-flight data. All of these comparisons indicate that the present 
method yields realistic results. In addition, the results obtained by the 
present method are considered to be more valid in most cases than the 
results obtained by the other methods and at least as valid in all cases. 
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INTRODUCTION 


A test conducted in a free-flight ‘ballistic range ordinarily yields 
the following basic data: 

1. The location of the model center of gravity in space as a 
function of time at a finite number of points 

2. The angles of attack, yaw, and roll as functions of time at 
a finite number of points 

From these basic data, other information such as the velocity as a function 
of time (or distance) and the period of the resulting oscillation can be 
determined . 

The problem of obtaining the aerodynamic coefficients from these 
data involves writing down the equations of motion and solving this set, 
if possible, using the given data. In the most general case, no solution 
to the equations of motion exists. It therefore becomes necessary to 
make certain assumptions, and hence simplify the equations to a set that 
can be solved. 

One assumption that is commonly made to simplify the equations of 
motion is that the restoring moment is linear with respect to angle of 
attack. Under this assumption, a number of solutions to the problem do 
exist (e.g., ref. l) . It has been found, however, that bodies of revolu- 
tion at high Mach numbers, especially blunt-nosed slender bodies, experi- 
ence extremely nonlinear restoring moments. It is thus necessary, in 
order to treat bodies of this type, to discard the assumption of a linear 
restoring moment and to investigate under what conditions and with what 
assumptions a solution to the nonlinear case can be obtained. 

A variety of solutions to this problem exist (e.g., refs. 2 and 3) • 
However, for a set of data obtained in the Ames Supersonic Free-Flight 
Wind Tunnel, the adequacy of these known methods was either doubtful or 
unknown. It was thus decided to attempt to gain better insight into the 
problem of analyzing the motion of a body governed by a nonlinear restor- 
ing moment. A systematic description of this investigation will be given 
in this report. Attention is first concentrated on a very simple nonlinear 
case. After examination of this solution, the development will proceed 
step by step to more complicated cases, eventually arriving at solutions 
that can be applied to free-flight data with a great deal of confidence. 


SYMBOLS 


A 


constant in equations (8) , ( 17 ) , and ( 25 ) 
constant in equation (l) 


B 


M3 4=r > 


3 


Cm 

c i%, 


restoring moment coefficient. 


M 

qS l 


moment-curve slope. 


dC m 

da 


Cma L 

F(cp,k) 

g 

h 

I 

K(k) 

k 

l 

M 

m 

n 

P 

<1 

R 

S 

T 

t 

u 


moment-curve slope obtained from linear theory (defined in 
eq. ( 6 )) 

incomplete elliptic integral of the first kind 

quantity introduced in obtaining the solutions to a number of 
integrals (defined in text for the various cases encountered) 

constant in equation (25) 

^moment of inertia about an axis through the center of gravity 
and normal to the axis of symmetry 

complete elliptic integral of the first kind 

modulus of elliptic integral (defined in text for the various 
cases encountered) 

reference length 

restoring moment 

constant in equations (8), (17) , and (25) 
constant in equation (17) 
exponent in equation (l) 

1 2 

dynamic pressure, — pV 
A 

qSZ 

reference area 
period of oscillation 
time 

factor appearing in expressions for co^ (defined in text for 
the various cases encountered) 


b 


V 

x 

Xi,X 2 

y 

yi>y 2 >y 3 

a 

a E 

a 0 

P 

P(a,t>) 

r(a) 

<p 

(•) 


velocity 



roots of the polynomial appearing in equation (19) 
a 

roots of the polynomial appearing in equation (26) 
angle of attack 

effective angle of attack (the true value of C at a = 
is equal to at a = am) 

maximum angle of attack 

(Note that when "reasonably planar motion” is considered^ 

% refers to the maximum resultant angle of attack.) 

minimum resultant angle of attack 

angle of sideslip 

beta function 

gamma function 

argument of elliptic integral (defined in text for the various 
cases encountered) 


dt 


( ) 


&t‘ 


( ) 


f) 
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SOLUTIONS TO THE PROBLEM 


Testing of bodies of revolution in the Ames Supersonic Free-Flight 
Wind Tunnel has indicated that the following conditions are satisfied in 
many cases: 

1. The roll rate of the model is small. 

2. The damping of the resulting motion is extremely small when 
a single cycle of oscillation is considered. 

3 . The resulting motion of the model is reasonably planar. (The 
definition that is herein adopted of "reasonably planar motion" 
is a motion which has the ratio o, 0 Ja^ <0.3 and which ro ughl y 
describes an ellipse in the a, 0 plane.) 

It will be assumed in the developments to follow that the roll rate and 
damping are identically zero and that the motion is identically planar. 
However, the solutions will be considered applicable if the conditions 
given above are satisfied. 


Moment Consisting of a Single Term in an Arbitrary Power of a 


A simple nonlinear case that we can consider is one where the moment 1 
can be represented as 


hence 


M = CmqS l = -Ba P p > 1 


( 1 ) 


C 


Bpa ?- 1 

“a “ “ qS l 


( 2 ) 


Given planar motion, with zero roll rate and damping, our equation of 
motion to be solved is 


la - M = 0 


or 


la + Ba p = 0 


It should be realized that the desired moment is an odd function 
[M( -a) = -M(a)] which does not follow from equation (l) . The derivation 
that follows, however, does not depend on this point. 
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On multiplying this equation by a and integrating, we obtain 

Id 2 Ba p+1 


2 p + 1 


= C 


To find C, we use the fact that a = 0 when a = c%. Thus 


C = 


Ba m J 


p+i 


p + 1 


and 


d = Jt^ti <^ +1 - ^ 


dt = 



da 

■K P+1 - cc . p+1 


( 3 ) 


We can determine T/k (where T is the period of the oscillation) by 
integrating equation ( 3 ) from a = 0 to a = Oj^; that is. 



l 


(p + 1)1 

2B 




da 


j 


p+i p+i 

a m - a^ 


This integration can he performed hy employing the beta function 


P(a,b) = J u a - 1 (l - up-Mu > 0 


Thus by letting 


nP+i 

u-r-sOi 


du = 


c^m J 

(p + l)a p (p + l)u p+1 da 




p+i 


da = 




we get 


T = / (p + 1) I r 1 

k v 2B 


Omdu 


0 , p?T / p+i p+i 

(p + 1) u \l oom - am u 


vo -<1 s=- > 


MJ -p- P> 
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or 


T = 


8l 


(p + ^Bctm 


FT ^ 


p + 1 ' 2 


We can now introduce the relationship "between the beta function and 
the gamma function that 


P(a,b) = 


T = 


r(a)r(b) 
r(a + b) 


8l 


+ 1 y(i 


(p + ljBom 


p-i 


P ± 3 . 

2p + 2 


and since 


r * = 


r(a + 1) = aT(a) 


our expression for the period reduces to 


T = 


2jtl(p + 3)' 


ZJt_2 


(p + 1 r MP + ^ 


Solving for B, we get 


B = 


2jtl(p + 3) 2 
(p + l)a m P T 2 






Putting equation (4) into equation (2) gives our final expression for 

2 


■'Hfet, 


2jtp(p + 3) : 

p + 1 


p + 2 


r ( 3£Ji_5 

V 2 P + 2 




p-i 




V T 2 qSZ 


( 5 ) 


and a similar expression can he obtained for C ; 


m* 
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For this simple case, our problem is thus solved. However, for the 
more complicated cases to follow that involve more than one unknown, we 
will have to develop an indirect method of finding the unknowns. Since 
this indirect method also leads to a better intuitive feel for the prob- 
lem, it will be developed for this case also. We proceed as follows: 
Assuming p = 1 (linear theory) in equation ( 5 ) , we get 


, _ -kiPl 

^ ' AST 


( 6 ) 


This is the familiar expression for Cmx under linear theory and our 
additional assumptions. Now, it can be argued that there exists an angle 
of attack (called a E ) through which our model has oscillated at which 
the local value of Cmx "the same as Cmx E obtained from linear theory. 

To obtain this angle of attack, we equate ( 5 ) and ( 6 ) , substituting a E 
for a in equation ( 5 ) • Thus 


— 2jrp(p + 3) 2 

p + 1 




P + 2 
P + 1 


3p + 5 

2P + 2 


2 



-4JT 2 ! 

I^qSZ 


from which 


r 

r ( 2 l ± 3 \ 

2 

2 jt(p + l) 

\ 2 p + y 


p(p + 3) 2 




L Vp + vj 

4 


(7) 


Now, equation ( 7 ) tells us that if our model was indeed governed by the 
restoring moment assumed in equation (l) (consider p as a known value), 
then we can determine the real angle of attack at which the Cmxp obtained 
from linear theory is equal to the local value of Cmx- In effect, the 
parameter ag/om tells us the transformation from a plot of Cmx E vs. am 
to a plot of Cmx vs. a. 


Figure 1 shows a plot of a-g/om vs. p. It is interesting to note 
that for a cubic moment (M = -Ba3) , CmxL> obtained using linear theory, 
should be applied at 0.489 <%• This is considerably lower than the root 
mean square value of a sine wave (O. 7 O 7 (%) which one might have expected 
to be a fairly close value. 

The downfall of the approximation that was made in equation (l) is 
that for p > 1 the moment-curve slope at a = 0° is 0. This is obviously 
too severe a restriction to place on our solution. However, it is felt 
that this single analysis has led to an intuitive feel for the problem at 
hand. We will now consider a slightly more complicated case. 


vo -fr > 


4 =-> 
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Moment Consisting of a Linear and a Cubic Term 


This problem has been treated by Rasmussen in reference 2. Rasmussen 
considered the two-degrees-of-freedom case (no restriction to planar 
motion) and obtained solutions in closed form. As a step in building up 
the present method of solution, however, it is necessary to treat this case 
in the same manner as our previous example. Also, certain approximations 
made by Rasmussen in simplifying his final results for application pur- 
poses can, for the case of planar motion, be applied in their exact form. 

It should be noted that the exact solutions presented in reference 2, when 
specialized to the case of planar motion, can be transformed through 
mathematical manipulations into the same solutions at which we will arrive. 

Consider the case where the restoring moment governing our model can 
be represented as 

M = Cm qSZ = -Aa - mAx 3 (8) 

hence 

Cma = “ il I C 1 + (9) 

This assumption is more realistic than our first example since, in this 
case, the moment-curve slope at a = 0° is not in general zero. It 
should be noted that under this assumption d 2 C m /dot , 2 = 0 at a = 0°. The 
equation of motion that will be solved is again 

Td - M = 0 
or 


Id + Aa + mAa 3 = 0 

Following exactly the procedure set forth in the first example, we obtain 



da 

+ monj 4 - 2a 2 - ma 4 ) 


f. 
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or on letting y = a/om 


T '= 


321 




w iJ ° /a 

\/ \ mam 


Am > 0 


> (10) 


T = 


321 


ay 


Am < 0 


-Amam J o ^ _ hk% 2 + 2 


mam 


Now we must consider the three possible cases that exist under equa- 
tion (8) . These are shown in sketch (a) . 





Sketch (a) 

Case 1.- Since Am is greater than zero, the first of equations (10) 
is the appropriate one. With m greater than zero, the polynomial under 
the radical within the integral in this equation is nowhere negative 
between the limits of 0 and 1, so we can proceed with the integration. 
This integration can be performed by using 213.OO from Byrd and Friedman 
(ref. 4) . 


T = 




mom 2 

2(mom 2 + l) 


where 


VO ->] > 


VO > 


n 


and K(k) denotes the complete elliptic integral of the first kind. 
Solving for A, we get 


A = g^k)] 2 (11) 

Putting equation (ll) into equation (9) gives our final expression for 




-l6[K(k)] 2 
mom 2 + 1 


(l + 3mx 2 ) 


I 

Q^qSZ 


( 12 ) 


Once again, we solve for the effective angle of attack at which the local 
value of Cn^ is the same as Cmxp obtained from linear theory. Hence, 
by equating equations (6) and (12), we obtain 


/ rt 2 ( 1 + m% 2 ) 1 

J 12m[K(k) ] 2 3m # 


( 13 ) 


Equation (13) tells us how to determine our effective angle of attack so 
that the value of obtained from linear theory is the local value 

of Cn^. In figure 2, ag/ccm is plotted as a function of raa m 2 . It is 
noted that the value of 0.489 obtained in the first example considered 
(where the linear term was neglected) is a limiting value in this case, 
and nowhere does the shift differ substantially from this value. 


Case 2 .- In this case, Am < 0, so the second of equations (10) is 
applicable. Before we can proceed with this integration, we must deter- 
mine under what conditions the polynomial under the radical within the 
integral is non-negative. If it is negative anywhere between the limits 
of 0 and 1, the period of the oscillation is imaginary, corresponding to 
instability. This polynomial can be shown to be non-negative between the 
given limits, leading to a period that is real and finite, as long as the 
following condition holds: 

-1 < m^ 2 < 0 (l4) 


This can be pictured more clearly be referring to sketch (b) 
tion = -1 gives the value of at 

which there is an unstable trim point; condition 
(l4) states that 0 < ocm < a^.. With these thoughts 
in mind, we can now integrate equation (10) by 
using 220.00 from reference 4. 


The condi- 


T = 


321 

-Amxni^ 


gK(k) 
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where 


g 2 = k 2 


-g % 2 
mo^ 2 + 2 


K(k) = complete elliptic integral 
of the first kind 


Solving for A, we get 


A = 


-321 

mom 2 ! 2 


Putting this equation into equation 


Cmx = -32[K(k )] 2 
mom 2 + 2 


g 2 [ K(k) ] 2 

( 9 ) gives the expression for 

(1 + 3mx 2 ) *- 

T 2 qSl 


Cma, 


The expression for ag Ts now obtained by equating this expression with 
equation ( 6 ) . 


/jt 2 (2 + mptm 2 ) 

•J 24m[K(k) ] 2 3m 


(15) 


This case is also plotted in figure 2. It is noted that ag/a m for this 
case lies between O. 5 O and 0.577* 


Case 3 .- In this case. Am > 0 so equations (10) assume the same form 
as under Case 1; that is, 


T = 



J (1 - 


d£ 



mom 2 + 
mctm 2 J 


Here the condition for the quantity under the radical within the integral 
to be positive or zero between the given limits can be written as 


or 


mot m 2 + 2 

mc%i 2 


> 0 


mom 2 < 


2 


(16) 


VO -^1 -r > 


VO -<] 4=r > 


13 


Physically, what this means can he explained as follows (see sketch (c)). 
If a model were launched at a = 0,4, it 
would continue to fly at that angle of 
attack throughout its trajectory. Simi- 
larly, if launched at a 3 , it would 
oscillate between <x 3 and 05 . The area 
under the curve between 0,4 and a 5 is 
a measure of the energy added to the sys- 
tem while the corresponding area between 
04 and a 3 is a measure of the energy 
removed from the system. Since we have 
neglected damping, the model would oscil- 
late between limits that make these areas 

equal. Now we have insisted that the model oscillate through 
(i.e., between a = and a = -a^) , which cannot occur unless 
large enough to allow the shaded area above -C m = 0 to be equal to that 
below -Cm = 0. 'Thus, condition (l 6) does nothing more than state that 
a m > a 2 . With this in mind, we find that the solution to Case 3 is 
identical to that of Case 1; that is. 



Sketch (c) 


o=0 
a m is 


I rt 2 (l + mom 2 ) 1" 

V 12m[K(k)] 2 3m 


(13) 


where 


k 2 _ mam 2 
2(mo m 2 + l) 

It is noted from figure 2 that for this case, 0^/0^ varies between 0.4o8 
and 0.489. As in Case 1, throughout most of the range the value obtained 
when the linear term is neglected (0.489) is a close approximation. 


USE OF LINEAR PLUS CUBIC APPROXIMATION 


Assume that at least two firings are made of a given configuration 
in a free-f light facility. Assume also that over one cycle of oscillation 
the damping is negligible, that the roll rate is small, and that the motion 
is reasgnably planar (see definition, page 5) . 

It is then an easy matter to determine and to apply linear 

theory and determine Crox L for each of the runs. If the data points 
appear to fall on curves like any of those shown in figure 3, the linear 
plus cubic approximation can satisfactorily fit the data. It is obvious 
that if the data show an inflection point or are other than continuously 
increasing or continuously decreasing, the linear plus cubic approximation 
cannot fit very well. 


Ill- 


Two runs are then chosen which have values of 0 % reasonably far 
apart. We* can then obtain two equations with two unknowns by inserting 
the following into equation (9) : Replace 

Cma *>y Cmxp 
a by a/g 

and define 
Then 

Cmx L ^ = -R(l + Bmo-E! 2 ) 

Cn&j^ = " R ^ 1 + 3ma.E2 2 ) 


In matrix form, these equations become 




1 30-E ! 2 



l 

I 

O 

1 

iT 

l 


1 3o-e 2 2 


mR 


or 


■CmxL 


P 


Q 


Since P is nonsingular, it has a unique inverse, and we can write 


1 1 

hd 

1 

H 

1 l 

_~ Cffla L 

= 

p-i 

P 

Q 

= 

Q 


We are thus led to the following expression for the unknown coefficients: 


R 


- a E s 2 ' 

O p 



a Ei - a-Ea 

mR 


1 


3( a El 2 - a E 2 2 ) 


a Ei 2 



2 „ 2 


- c ma Ll 

a Ei " a E 2 



-1 


-Pm 

3(a-E x 2 " a E 2 2 ) 


^Ie 


It will be obvious from a plot of Cmxp vs * “m whether the test results 
qualify for Case 2 or for Cases 1 and 3 (see fig. 3 ) . 


VO -^1 4^ > 


< -*■ r- o\ 
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After choosing the governing case, it then becomes merely a matter 
of assuming a value of a E (anything near 0.5 a m is sufficiently close), 
determining the P" 1 matrix using equation ( 13 ) or ( 15 ) , and from this 
forming the Q matrix. The quantity m is obviously mR/R. 

This new value of m is used to obtain a new P -1 matrix and the 
iteration process continues. The iteration process converges very rapidly 
for each of the three cases. After m and R are obtained, C m , as a 
function of a, is determined from equation (8) . 

After the iteration process has converged, the solution will of 
necessity go through the two data points used. The demonstration of the 
adequacy of the linear plus cubic moment assumption is in how well the 
curve goes through additional data points obtained from testing the same 
configuration. This is checked by constructing the Cn^ vs. a m curve 

from the following equation, using the values of R and m given by the 
iteration process: 


- “R(l + 


If this leads to a reasonable fit of the data points, a better fit can be 
obtained in the following manner. First a least squares fit is obtained 
of Cn^ vs * 0/111 ^ a ' v 'i n S the same form as equation ( 9 ) • 

T 

= a + bOm 2 

After solving for a and b, two points are chosen from this least squares 
fit to use as inputs for the iteration procedure. The iteration procedure 
will then yield a fit of the data points that falls very close to the 
least squares fit. This could best be termed a quasi least squares 
procedure. 


Moment Consisting of a Linear, a Cubic, and a Quint ic Term 


Next, consider the case where the restoring moment governing our 
model can be represented as 


M = CmqSZ = -Act, - mAx3 - nAa 5 


(17) 


l6 


hence 


CiAx 


^ (1 + 3m. 2 + 5 m‘) 


(18) 


Once again it is noted that 
that will he solved is 


— = 0 at a = 0° . 

da 2 


The equation of motion 


la - M = 0 


or 


la + Aa + mAa 3 + nAa 5 = 0 

Again we follow exactly the procedure set forth in the first example 
considered and obtain 


T = 4~9 61 




da 


sj A (6am 3 + 3 mam + 2nam 6 — 6a 2 - 3^^ 2na 6 ) 


or on letting x 



T = 


121 
am 4 v 


dx 


An x( 1 - x) 


+ 


3 mam + 2na m " 
2nam 4 


x + 


6 + 3ma ia 2 + 2na m " 
2nam 4 


']} 


( 19 ) 

We will define by x x and X 2 the roots to the quadratic polynomial 
that appears in the denominator of our integral. 

-(3m + 2nam. 2 ) + J 9m 2 - W 3 n - 12mnam 2 - lSin 2 ^ 4 

Xi = 

4na m 2 

-(3m + 2m ra 2 ) - J~9ra?- - 46n - lSnnam 2 - 1 2n 2 a m 4 

X2 = 

4nam 2 


These roots may be both real or both complex. Now, equation (19) may be 
rewritten in the following manner. 


vo -j -fr- > 


17 


T = 


121 


dx 


Amm 7 x(l - x) (x - xi) (x - x 2 ) 


An > 0 


I 121 ^ dx 

7 -Ana® 4 J 0 ,y_ x (i - x )(x - x x )(x - x 2 ) 


> ( 20 ) 

An < 0 


The quantity under the radical outside the integral is then always posi- 
tive. To obtain a meaningful solution to equations (20) , it is necessary 
for the quantity under the other radical to he non-negative between the 
given limits. This condition will be discussed when the individual cases 
that exist under equation ( 17 ) are considered. These cases are shown in 
figure 4. 


Case 1:^ Linear, cubic, and quintic terms are a ll stabilizing (A > 0, 
m > 0, n > 0) .- The first thing we must determine Is nndp-r -irba-h 
the polynomial under the radical within the integral in the first of equa- 
tions (20) is non-negative between the limits 0 anri 1 . 

Assume first that x ± and x 2 are real (i.e., 9m 2 - 48n - 12mna m 2 
- 12n 2 a m 4 > 0) . This can obviously be achieved by making m large and 
n small. Under this assumption, it is apparent that x x > x 2 and that 
x 2 < 0. To show that xi < 0, rewrite the expression for xx as follows: 

_ - 7 9^ 2 +12mnct m 2 +4n 2 o m 4 + J ( 9m 2 +12mna m 2 +4n 2 a m 4 ) - ( 48n+24mna m 2 +l6n 2 a m 4 ) 

Ana m 2 

Since the expression under the first radical is larger than that under 
the second radical and the denominator is positive, xi < 0. Wow, since 
x 2 < x i <0.? the terms (x - Xi) and (x - x 2 ) are necessarily positive 
between the given limits and the expression under the radical is >0 
between these limits. We thus have the following (see sketch (d))7 



Sketch (d) 


l8 


Next,, assume that x x and X2 are complex ( i .e . , 9m 2 - - Wn - 12mncc m 2 
- 12n 2 a m 4 < 0) . This can obviously be achieved by making n large and 
m small. 

Now, Xi and X2 can be expressed as follows: 

-a + i \fb 
x x = 5 

-a - i >v/b 


with a > 0, b > 0, c > 0. Then the product (x - x x ) (x - X2) has the 
following expansion: 


/ w \ p 2ax a 2 b 

(x - x x ) (x - x 2 ) = x 2 + -5- + + ^2 

This term is positive for 0 < x < 1, and hence the term in the denominator 
of our integral is >0 between the given limits. We thus have the fol- 
lowing (see sketch (e)): 



We can now proceed to solve equations (20) . The two subcases of 
Case 1 will be treated separately. 

Subcase 1, x x and x 2 real: Making use of 256. 00 from reference 4, 

we can write 


T = gK(k) 

n/ Ana^ 


19 


where 


g 2 = 


4 

XxX 2 - X 2 


k 2 = X 1 ~ x 2 
’ X X X 2 - x 2 


Solving for A, we get 


A = 


121 

no m 4 ! 2 


g 2 [K(k)] 2 


Putting this equation into equation (l8) gives the expression for Cop, 


Cap, = - 


12g 2 [K(k) ] 2 
“% 4 


(l + 3ma 2 + 5 r x<< 4 ) 


I 

l^qSZ 


The expression for Oj. is now obtained by equating this expression with 
equation (6) . 


a E 


/-3m + 79m 2 - 20n(l - u) 
7 lOn 


( 21 ) 


where 


a 2 rK % 4 

3lf[ K(k)] 2 

The plus sign in front of the inner radical in equation (21) is chosen in 
order for co^ to be real. 

Subcase 2., Xi and x 2 complex: Applying 259*00 from reference 4, we 

can write 


where 


T = 


I 48l 
Anc% 4 


gK(k) 


7 XiX 2 (l - Xi - X 2 + XiX 2 ) 

mm 4 

J \ (1 + mom 2 + nom 4 ) (6 + 3mOm 2 + 2no m 4 ) 
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k 2 = 




Xi + X2 - 2XiX2 

4 ^X 1X2(1 - Xi - X2 + X1X2) 


1 2 

2 " 4 


(4 + 3ma^ + 2mO_ 


■J 6(l + mam 2 + nc^ 4 ) (6 + 3nr%i 2 + 2nc% 4 ) 


Solving for A, ve get 


A = 


48 1 


g 2 [K(k) r 


□am 4 ! 12 

Fatting this expression into equation (18) gives the expression for Cn^. 


C nbc, “ 


-48g 2 [K(k) ] 2 


(l + 3rax 2 + 5na 4 ) 


nom** T*q.Sl 

Equating this expression with equation (6) gives the expression for ag. 


a- 


E 


I -3m + ^9m 2 - 20n(l - u) 

__ 


( 22 ) 


where 


j^nom 4 

12g 2 [K(k) ] 2 

It is noted that equation (22) has the same form as equation (21 ) , hut it 
should he kept in mind that u, g, and k are expressed differently in the 
two subcases. 

General discussion concerning additional cases .- A rigorous deriva- 
tion of the expression for ot/g for each of the additional cases that 
exist under equation (17) would consume much space and serve little 
purpose. The form of this expression for each case can he written as 


a E 


/ -3m ± s/9m 2 - 20n(l - u) 
s lOn 


(23) 


The appropriate expressions for u, g, and k vary depending on whether 
xi and X2 are real or complex and also on the case under consideration. 
We will define a number of u, g, and k which will allow us to. indicate 
in tabular form the appropriate expressions to use. 


VO -O -p- > 
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Case 


1 


2.1 


2.2 


3 


Ui = 


gi 2 = 


SgStKCk )] 2 

4 

I XiX 2 - x 2 I 


U2 = 


-j^noj,! 4 

3g 2 [K(k) f 


S 3 2 = 


nc%n 


u 3 = 


S 2 2 = 


n 2 mm 4 

12g 2 [K(k) ] 2 

4 

|XiX 2 - Xil 


y| ( 1 + rattm 2 + ntXjjj ) ( 6 + 3^^111 ^nctjjj ) 


ki = 


X 2 - X! 
XiX 2 - x 2 


k 2 = 


x 2 - Xi 

XlX 2 - Xi 



3(4 + 3nK% 2 + ^Om 4 ) 

kj 6(l + mom 2 + nam 4 ) (6 + 3m m 2 + 2na m 4 ) 


TABLE I 


Description 




Valid 

root 


xi and x 2 

u 

s 2 

k_ 

(eq. (23)) 


a li a 2 

Real 

ui 

Sl 2 

ki 





or 




+ 




Complex 

U3 

S3 2 

k 3 





Complex 

only 

U 3 

g 3 2 

k 3 

+ Always 
- Sometimes 




Real 

or 

Ul 

g! 2 

ki 

+ Always 
- Sometimes 

O-l = , 

1- 

- Jm 2 -bn 
2n 

Complex 


P 

k 3 

a 2 = . 

l-m 

£ 

i 

S' 

CVJ 

+ 

^3 

g 3 


4 

2n 

Real 

only 

U2 

Si 2 

ki 

- Always 
+ Sometimes 

*1 = , 

/- 

- Vm 2 -4n 

2n 


a 1' a 
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TABLE I - Concluded 



Ti 


A<0 

m<0 

n<0 


Valid 

root 

X! and x 2 _u_ _k_ (eg. (23)) 


a X^2 


Real g x 2 k x 

only 


Real ui g 2 2 k2 

or 

Complex U3 g 3 2 k 3 


a i 


a 2 


-F 1 


Jw£-kn 


2n 


/ -3m - sj9m 2 -kQn 
J kn 



Real m g 2 2 k 2 

or 

Complex u 3 g 3 2 k 3 


a 2 


■F 


3m - \/9m 2 -48n 


4n 


•C« 


'aTa\ 


A<0 
m< 0 
n>0 


2 a |\ a 


Real U2 g 2 k2 Both 
only always 


a i 


a 2 




\/m 2 -An 


2n 


/ -3m - 

V An 


The column headed "valid root" denotes the appropriate sign to he 
used in front of the inner radical in equation (23) . Some cases can have 

two valid roots because of the double- 
valued nature of Cu^. It is clear 
from sketch (f) that there are two 
transformations that can he applied 
to point A that will place it on the 
Cma vs. a curve, while only one such 
transformation exists for point C. 

The point labeled B is evidently the 
dividing point between these two 
possibilities . 



The tabulated values of and 
a 2 indicate the dividing line between 
Sketch (f) regions where our solution does and 

does not exist. Specifically, we must 
insist that <x m < a,i and/or a m > a 2 . If these conditions are not satis- 
fied, the quantity under the radical within the integral in equations (20) 
becomes negative and the solution that has been developed has no meaning. 
The physical meaning of these conditions can be realized by noting that 


<$ -=i* ON 
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we have insisted that the model oscillate through a = 0° and, if the 
conditions are not satisfied, this will not he the case. It should he 
noted that a x is an unstable trim point and that a 2 is located so 
that the shaded areas above and below -C m = 0 are equal. 

USE OF THE LINEAR PLUS CUBIC PLUS QUINTIC APPROXIMATION 


Before this approximation can be applied, at least three firings of 
a given configuration must have been made. We then proceed as in the case 
of the linear plus cubic approximation (page 13 ) and obtain three equa- 
tions with three unknowns. These equations can be written in matrix form 
as follows: 


-Cmcx,^ 


1 3<* El 2 5 o- E i 4 


R 

- Cm ^L 2 


1 3^e 2 2 ^Es 4 


mR 

" c ^L 3 


1 3o-e 3 2 5°-E 3 4 


nR 


Once again we determine the inverse of the 
then write 


square matrix, and we can 


R 


mR 


nR 


a E 2 2a E 3 2 

(a El 2 - a Es 2 ) (a El 2 - a Ea 2 ) 

1 a E 2 2 + S'Ea 2 

3 ( a Ei 2 “ a E 2 2 ) (* El 2 " a E 3 2 ) 


5(a-E:L 2 - a E 2 2 ) ( a E x 2 - a E 3 2 ) 


-^i^s 2 

(a El 2 - a E2 2 )(a E2 2 - a Es 2 ) 

1 ^ 

3 ( a Ei 2 " a E 2 2 ) ( a E 2 2 “ a E 3 2 ) 


5(a E i 2 - a Es 2 ) (a E2 2 - a E3 2 ) 


^i 20 ^ 2 


- C »aL 1 

(aEi 2 ” a E 3 2 ) ( a E 2 2 - a E 3 2 ) 


1 ' W * ^ 


.C- 

3 ( a Ei 2 " a E 3 2 ) ( a E 2 2 " a E 3 2 ) 



1 


- c ”-l 3 

5(a.E 1 2 - a E 3 2 ) ( a E 2 2 “ a E 3 2 ) 



(24) 


From equations (2k), we would like to determine R, m, and n. 

From a plot of Cn^ vs. (% and reference to figure k, we can 

L 

determine the cases under which our test results may fall. The solutions 
are so similar that it need not concern us that we are unable to pinpoint 
a specific case that governs our data. This is especially true if the 
problem is programmed for a digital computer. It is fairly easy to set 
up a general program that will handle all possible cases. 
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A guess is now made at values of m, n, and R (only a positive or 
negative sign need be assigned to R) and an iteration process is started. 

From m ; n, and R, the three 
values of cc/g can be determined 
a m from equation (23) and the infor- 
mation in table I. When more 
than one root of equation (23) 
is valid , the root that gives 
points -the largest spread between the 

highest and lowest values of a,g 
should be used. Sketch (g) clar- 
ifies this. It is apparent that 
the points A, B, and C better 
define the Cmx vs. a curve than 
do the points A , B, and D. Equa- 
R , mR , and nR, and hence for m and n. 
are used to obtain new values of the 



O— Doto 


Sketch (g) 

tions (24) can then be solved for 
These new values of m, n, and R 


cc-e, and the iteration process continues. 


It has been found that the iteration process converges very rapidly 
for all cases that are initially stable (Cases 1 through 4) , even with 
extremely poor initial guesses of m and n. The same is unfortunately not 
true for cases that are initially unstable (Cases 5 through 7) . For these 
cases the iteration process diverges, even with almost exact guesses of 
m and n. This condition can be remedied to some extent as follows: We 

will make the assumption that R (i.e., -Cntc, at a = 0°) is a known value 
and will then proceed to obtain expressions for m and n. Only two data 
points will be used to get these values; the third data point will be used 
in a way that is explained presently. 


Regarding R as known, our two equations are 


R 

^L 2 


+ 1 


+ 1 


3 a-Ei 2 

5 %; 

CM 

& 

50 ^ 2 ' 


m 


n 


By determining the inverse of the square matrix that appears above, we 
can write 


m 


a E 2 2 “O-Ei 2 


" P *Wa . " 


3 a E 1 2 ( a E 1 2 “ a E 2 2 ) 3 a E 2 2 ( a E 1 2 “ a E 2 2 ) 


+ x 

■n 


-1 1 



11 


_5a El 2 ( a 'E 1 2 - a E 2 2 ) 5a-E 2 2 ( a E 1 2 - a E 2 2 )_ 


R + 1 


vo -j -t=- > 


M3 -J > 
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Now we can make a guess of m and n and go "through an iteration process 
as before. In this case the convergence is rapid if the assumed value 
of R is anywhere reasonably close to its true value. 

This process is repeated for several different assumed values of R. 
In each case, after the iteration has converged, we form the following 
parameter: ACmxp = CnttL (known 

for third data point) - Cn^, 

( given by iterated solution). 

Then we can graphically deter- 
mine a value of R so that the 
iterated solution will go 
through all three data points, 
and with this value of R as 
an input, the iteration will 
yield the corresponding values 
of m and n (see sketch (h)). 

Regardless of which case is under consideration, once the iteration 
process has converged, the solution will of necessity pass through the 
three given data points. The adequacy of the 1-3-5 approximation is deter- 
mined by seeing how well it fits additional data points obtained from 
tests of the same configuration. The Cmxp vs. a m curve is constructed 
from the following equation using the values of R, m, and n given by 
the iteration process: 



Cmx, L = -R(l + 3maE 2 + 5na-E 4 ) 

Once again, a better fit can be obtained by using the quasi least squares 
procedure mentioned on page 15 • For this approximation, a least squares 
fit of the form 


Cjoxl = a + t^m 2 + com 4 

is obtained, and then three points from this fit are used as inputs for 
the iteration procedure . 
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Moment Consisting of a Linear, a Quadratic, and a Cubic Term 


Next consider the case where the restoring moment 2 governing our 
model can he represented as 


2 3 

M = CmqS l = -Aa - hAa - mAa 


(25) 


hence 

= itl (l + 2ha + 3ma 2 ) 

For a moment of this type, it is noted that d 2 Cm/da 2 is not, in general, 
zero at a = 0°. The equation of motion that will he solved is 


Icg + Ax + hAx + mAa = 0 


Again the procedure set forth in the first example that was considered is 
followed, and we obtain 


T = -J 961 


'am 


da 


o J A(6a m 2 + 4hXm 3 + 3 mam 4 - 6a 2 - Ij-ba 3 - 3ma 4 ) 


or on letting y = a/am 


32i g f 1 dy 

Araxm '■'o IT » [" o / 3 mm 2 + 4ham \ Z f 3mxm 2 + 4ha.m + 6\ Amam 2 + ^ham + 6V | 

J ^ L \ 3raain 2 ) \ 3mam 2 ) ^ \ 3mam 2 /J 


T = rj&. r 1 





«/ 


(26) 


2 It should be realized that the desired moment is an odd function and 
that to be strictly correct we should write a| a | for the quadratic term. 

This does not make any difference in our derivation and hence will not be 

carried through. 


vo -j 4=- > 


M3 -J fr > 
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We will define "by y x , y 2 , and y 3 the roots of the cubic polynomial that 
appear in the denominator of our integral. Since these roots may all he 
real or one may he real and two complex, they will he defined as follows: 

# 

If there is only one real root, this root is defined as y x . 

Wo distinction between y 2 and y 3 is necessary . If there are 
three real roots, they are defined such that y x > y 2 > y 3 . 

To obtain expressions for the three roots, it is convenient to define the 
following quantities: 


4h 

c = — + 1 


d = 


3mcqja 
6 

3mam^ 

1 3 _ 2 


G 0 = ^ (2c 3 - 9c + 27 c - 9cd + 27 d) 


Ho = ^ (3c + 3d - c2) 


G = 3 J -G 0 + Jg 0 2 + Bo 3 
H = 3/ _G 0 - / G 0 2 + Ho 3 


(I = cos -1 ( + 


' -Gr 


Ho 


- if 

Go 

> 0 

+ if 

Go 

< 0 


0° < p. < l 80 ° 


Then, if Go 2 + Ho 3 >0, there are one real root and two complex roots 


yjL = G + H - J 


1 

y 2,3 - 2 


-(G + H + y) ±i ^3(G - H) 


(27) 


2. 3 

If Gq + Ho <0, there are three real roots 


/ \ i = 1 when p = 0 

Yi = “ v + 2 sf - Ho cos (y + 120°p J i = 2 when p = 2 

^ ^ i = 3 when p = 1 


(28) 
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Nov by making vise of 253*00, 255*00, 257*00, 259*00, and 260.00 from 
reference lj-, ve can obtain solutions to equations (26). The possible 
cases that exist under the 1-2-3 approximation are the same as under the 
1-3-5 approximation (fig* 4), with the exception that d 2 C m/da 2 is not 
generally zero at a = 0°. The solution for the effective angle of attack 
for each of the cases of interest can be written in the following form: 


-h ± s/h 2 - 3m(l - u) 
a E 3^ 


(29) 


With the aid of the following definitions, the appropriate expressions 
for u, g, k, and 9 to be used for a particular case can be presented in 
tabular form. 


Define 
gi 2 = • 


ga 


S3 


g4 


1 

J ( l^y 2 -y3+y 2 y 3 ) ( y^-yiYs-y 

t ' 

(yx ~ i) (y 2 - y 3 ) 

4 I 

(y 3 - T)(y2 - yi) | 

k 

(y2 - i)(y 3 - yi) 


yi^ 2 , 3 from equations (27) 


y x 2 3 from equation (28) 


k, 2 = l - 


2yi-y 2 -y 3 -yiy 2 -yiy 3 +2 y2y3 

4 J ( i-y 2 -y3+y 2 y3) ( yi^yiyg-yiys+ygya) 


k2 2 = 5 + 


2yi-yg-y3-yiy2-yiy3+2y 2 y 3 

4 7(i- y 2 -y 3 +y 2 y3) (y^-yiys-yiys+ygys) 


) y 1)2 3 from equations (27) 


2 _ (y 3 - i)(y 2 - yi) 

(yx - 1) (y 2 - y 3 ) 

2 = (y 2 - i)(y 3 - yi) 

4 " (y 3 - D(y 2 - yx) 

2 _ (yx - i)(y 3 - yg) 

5 " (y 2 - 1) (y 3 - yx) J 


\ 

from equation (28) 


yx,2, 3 


VO -<l -1 ~ > 


vo -a -t=- > 
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<Pi = 
cp 2 = 

9s 

cp 4 

9s 

ui 

^3 

Case 


2.1 


I J yi 2 -yiy 2 -yiy3+y2y 3 + yi •A-y 2 -y 3 +y 2 y 3 \ 

S x I ■ 1 1 I 

VVy^-yiya-yiya+yaya - yi Ji-yz-y^y&aJ 
s _i /-V yi 2 -yiy 2 -yiy 3 +y 2 y 3 + y x >/i-y 2 -y 3 +y 2 y 3 \ 

\ 7yi 2 -y lya-yiya+ysya + yi J^-yl-y^y^ J 


» y x 2 3 from equations ( 27 ) 


= sm 


- 1 / 


y3 - y2 

y 2 (ys - 1) 


= sm 


- sm 


-1 


! y 2 - yi 

yi(y 2 - 1 ) 


yi - ys 

ys(yi - 1 ) 


y 1 ^ 2 ^ 3 from equation ( 28 ) 




8g 2 [2K(k) - F(<p,k)]‘ 

8g 2 [F(9.k)] 2 . 

rn^mam 2 

8g 2 [F(cp,k)] 2 


Description 


-C n 


f A>0 
h >0 
m>0 


’Cm 


TABLE II 


yx>Jz>yz 

(i) 


yi < o 

y 2 ,y 3 complex 


Valid 

root 

_u_ _gf _k_ (eg. (29)) 

ui gi 2 ki cp^ 


y 3 < y 2 < yi < 0 ua g 4 t 5 <p 5 


A>0 . / 

h<0 / / yi < 0 

"^9/ /^Limiting ^,y 3 complex 
^~v^/ co se 


Ui gi 2 ki <Pi 


+ Always 
- Sometimes 
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TABLE II - Concluded 


Case Description 


2.2 


"Cm 


"Cr 


y±>y2>ys 

(i) 


yi < o 

y 2 ,y 3 complex 


Valid 

root 

JL Jl JL (eg. (29)) 


ai ,a 2 


Si 2 ki (Pi 


- h-^h 5 ^ 

+ Always a i - ^ 

- Sometimes 


a^a 2 a 

yi > y 2 > i 

y 3 < 0 

U3 

g 2 2 

^3 


A>0 

yi > 1 

y 2j ,y 3 complex 

u 2 

Si 2 

k 2 

<P 2 

/^\ h >0 

y Wo 

or 





a) d 

yi > 1 
y 3 < y 2 < 0 

ua 

S3 2 

k 4 

9 4 

A>o 






v h<o 
\m<0 

yi > 1 

y3 < y 2 < 0 

Us 

g 3 2 

k 4 

<P 4 


- Always 
+ Sometimes 


cc,! = 


-h- 7 h 2 -4m 


2m 


a x 


-h- V h 2 -4m 


2m 



yi < 0 

y 2 .y 3 complex 


yi < o 

y 2 ,y 3 complex 


yi > i 

y 2 ,y 3 complex 


^1 gl 2 ki cp! 


gl 2 ki CP! 


2 , Both 

^ Sl *2 <P 2 always 


a2 = 


-2h- 7^h 2 -l8m 
3m 


a 2 = 


-2h ~ 4 kh 2 -lQm 
3m 


a x = 


a 2 - 


-h+ \lh 2 -bm 
2m 

-2h- s/ 4h 2 -l8m 

3m 


(i) 


Everything indicated in this column is correct, hut the indicated location of the roots is 
not necessarily complete. However, all possibilities concerning the roots are covered 
(except yi > y 2 > y 3 > 1 which cannot occur) so an additional possibility under any 
given case can be analyzed by referring to one of the other cases. 


( 2 ) 


a 2 = 


_ -h + 3/h 2 - 4m + 2 s/h 2 - 3hVh 2 - = 


4m 


6m 


This table has been set up in the same manner as table I (the 1-3-5 
polynomial approximation) and a discussion of the significance of the 
various columns will be found immediately following that table. 


VO -vl 4=- > 


VO .f=- > 
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USE OP THE LINEAR PLUS QUADRATIC PLUS CUBIC APPROXIMATION 


The use of this approximation is very similar to use of the 1-3-5 
approximation previously discussed, and thus only the pertinent equations 
will he presented. 


•D 



a E 2 a E 3 

■" a Ej a 'E 3 

®Ei®EJ2 ” 


_r< 

X\ 


( a Ei - 

• a-E 2 ) ( a E x " a E 3 ) 

(a El - a Ea ) (ag 2 - a Es ) 

(a El - ag 3 ) (a Ea - a Ea ) 


tmx Ll 

hR 



'( a E 2 + a E 3 ) 

a Ei + ^ 

+ a E 2 ) 


~ Ca a-L z 


2(a El 

- <*E 2 ) ( a E! - a E 3 ) 

2(a El - a Es ) (ag 2 - a Es ) 

2(a El - a Es ) (a E2 - a Es ) 


mR 



1 

-1 

1 


- Cn b-L 3 


3( a 'E 1 

- a E2 ) (a El - 

3 ( a E! “ a E 2 ) ( a E 2 " a E 3 ) 

3(a El - a E3 ) (a Es - a E3 ) 



where ag is obtained from equation ( 29 ) . 

Once again an iteration process is used to obtain values of R, h, 
and m. This program has been tried for configurations that were initially 
stable (Cases 1-1)-) and was found to converge rapidly. It has not been 
tried for configurations that were initially unstable, but it is assumed 
that the iteration process will diverge as it did for the 1-3-5 approxi- 
mation. If this is the case, the following equations and the scheme 
indicated in sketch (h) should lead to convergence. 


h 


a E 2 

-a El 


c « Li 



2a Ei( a E 1 - a E2) 

2a E 2 ( a E 1 “ a E 2 ) 


R 

m 


-1 

1 



in 


3 a Ei( a Ei - a E 2 ) 

3 a ’E 2 ( a Ei - a E 2 ) 


R +1 


RESULTS AND DISCUSSION 


Tests conducted in the Ames Supersonic Free-Flight Wind Tunnel of a 
variety of different models have yielded the following interesting infor- 
mation concerning linear theory. It has been experimentally determined 
that the application of linear theory to a given model at a given maximum 
resultant angle of attack will, even for models governed by nonlinear 
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moments, yield a value of Cma,L "that, with certain qualifications, does 
not depend on the type of motion encountered. These qualifications are 
that the ratio of coo/oc.m should "be less than 0.3 (there is no sudden 
transition when oc,o/o«m = 0.3, "but above this value the preceding statement 
becomes questionable) , and the resulting motion in the a, 0 plane should 
roughly describe an ellipse. When these two conditions are satisfied, 
the motion is (arbitrarily) defined as reasonably planar motion, and the 
solutions developed in this report (given in addition small roll rate and 
damping) will apply. An example of reasonably planar motion is shown in 
sketch (i) and examples of nonplanar motion are shown in sketches (j) 
and (k) . 



A number of shots of a particular blunt-nosed body of revolution 
were made in the Ames Supersonic Free-Flight Wind Tunnel. The data points 
obtained from these shots by applying linear theory are shown in figure 5. 
The ordinate is the moment-curve slope from linear theory and the abscissa 
is the maximum resultant angle of attack. Each data point represents one 
shot of the given configuration and all data points were obtained at a 
Mach number of approximately 11. For each of these shots the ratio of 
a o Awn was less than 0.2, and for the majority of the shots this ratio was 
less than 0.1. It was hence assumed that the condition of reasonably 
planar motion was satisfied. Since the roll rate was small and the damp- 
ing practically zero over one cycle of oscillation, it was concluded that 
all of the foregoing conditions were satisfied, and the results of the 
derivations presented in this report could be applied to this collection 
of data. 

The first thing noted from figure 5 is that the moment-curve slope 
obtained from linear theory appears to get smaller with increasing maximum 
resultant angle of attack up to about 9° and then increases drastically. 
This type of behavior eliminates the 1-3 approximation from consideration, 
since the 1-3 approximation is limited to either a continuously increasing 
or a continuously decreasing moment-curve slope. However, it is apparent 
from figure 4 that our data appear to fall under Case 2.1 and that both 
the 1-3-5 approximation and the 1-2-3 approximation have a chance of 
adequately fitting the data. 
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The iteration type of solution that has been described was used to 
obtain both a 1-3-5 approximation and a 1-2-3 approximation to the given 
data. These approximations are shown in figure 6. The quasi least squares 
procedure mentioned on page 25 was used to obtain the fits. It is clear 
from figure 6 that the 1-3-5 approximation fits the data very well. It 
is also evident that it is much superior to the 1-2-3 approximation in 
fitting the data (the two fits had least square values of O.OOO985 and 
0.00509, respectively). It is felt that this latter conclusion will apply 
to the majority of configurations encountered since the behavior of most 
configurations at low angles of attack is fairly linear. 

Figure 7 presents curves of the moment coefficient plotted against 
the true angle of attack corresponding to the two approximations shown in 
figure 6. Although it has been concluded from figure 6 that the 1-3-5 
approximation is superior to the 1-2-3 approximation, it is interesting 
to note that the two curves in figure 7 show the same general trends and 
throughout most of the angle-of -attack range agree fairly closely with 
each other. 

In figure 8, a comparison is made between free-f light data reduced 
by using the 1-3-5 approximation and data obtained from tests conducted 
in the AEDC B-Minor Wind Tunnel. Curves are shown of the moment coeffi- 
cient plotted against angle of attack at Mach numbers of 5 and 11 from 
free-flight tests of a given configuration. A corresponding curve at Mach 
number 8 obtained from wind-tunnel tests is also shown. If it is assumed 
that the variation of the moment coefficient with Mach number at a given 
angle of attack is fairly linear, then the wind-tunnel data would be 
expected to fall about halfway between the two curves of free-flight data. 
This is essentially borne out throughout most of the angle-of-attack range. 
It should be noted that the configuration tested in the wind tunnel was 
very similar, but not identical, to the configuration tested in the free- 
flight facility. A small difference in the nose shape between these two 
configurations would be expected to modify the moment coefficient. It is 
thus felt that the method of reducing the free-flight data by the 1-3-5 
approximation has led to a realistic moment curve. 

In figure 9, the moment coefficient as a function of angle of attack 
is shown as obtained from a given set of free-flight data by using the 
1-3-5 approximation and by using the method of Rasmussen (ref. 2) . The 
curve showing the 1-3-5 approximation is reproduced from figure 7* The 
other curve was obtained as follows: A plot was made of vs. 

(% 2 + a 0 2 ) . For the method of reference 2 to be directly applicable, 
this plot should result in a straight line. This was not realized for 
the given set of data, but is was possible to draw three straight lines 
that came close to passing through all of the data points (the three lines 
were determined by least squares fits) . The method of reference 2 was 
then applied individually to each of the three line segments, and the 
three corresponding segments shown in figure 9 were obtained. It is noted 
in figure 9 that the two methods yield values of the moment coefficient 
that remain fairly close to each other throughout the angle-of-attack 
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range. It is felt, however, that the 1-3-5 approximation is the more 
applicable of the two methods to the particular set of data considered 
because of the nonlinearity of the Cn^ vs, (am 2 + ao 2 ) plot and the 
accompanying segmented approximation. 

Results obtained by using the 1-3-5 approximation are next compared 
with the corresponding results obtained by using the method of Murphy 
(ref. 3) . To apply this method, the restoring moment is assumed to be 
linear and the motion in the a, p plane for each run is approximated 
by an epicycle with damping (two rotating vectors, the tail of one rest- 
ing on the head of the other) . An arbitrary nonlinear moment approximation 
in odd powers of the resultant angle of attack can then be assumed 
(i.e., M = ki<x + k^a 3 + . . .) and solutions for the ki are obtained 
from knowledge of the frequencies and lengths of the rotating vectors 
found in the linearized solution. For a given run, the frequencies are 
constants and the lengths are approximated by their values at the midpoint 
of the trajectory. 

Several points about this method should be discussed before 
proceeding: 

1. The solutions obtained by Murphy are not exact from a mathematical 
viewpoint, except for the case of a linear restoring moment. The assump- 
tion that an epicycle is the solution to the equation of motion, when the 
restoring moment is nonlinear, leads to terms which contain mixed frequen- 
cies, and these terms are neglected. The resulting solutions can best be 
described as being the first step in an iteration process. 

2. When the 1-3 solution obtained by Murphy is applied to the 
corresponding case of planar motion, excellent agreement with the exact 
solution is obtained in most cases. In fact, Murphy’s solution corresponds 
to the line a^/am = 0.5 in figure 2, which would lead to reasonable 
results except in the regions -1 < mam 2 < -0.8 and -5 < mam 2 < -2. In 
these regions the error builds up very rapidly, reaching 100 percent in 
predicting the period of the resulting oscillation (finite versus infinite) 
at the points mam 2 = -1* -2. 

3* When the 1-3-5 solution obtained by Murphy is applied to the 
corresponding case of planar motion, good agreement with the exact solu- 
tion is obtained in many cases. Because of the number of variables 
involved, however, it would be very difficult to specify the range of 
applicability of this approximation as was possible for the 1-3 case. 
Comparisons have indicated that for this case it is not uncommon for 
Murphy 1 s method to differ from the exact solution by 10 percent in 
predicting the period of the oscillation. 

In figure 10, the moment coefficient as a function of angle of attack 
is shown as obtained from the given set of free-flight data by using the 
1-3-5 approximation of the present report and by using the method of Murphy 
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(ref. 3 ) . The 1-3 approximation of reference 3 has been used and applied 
in segments, each segment being determined by a least squares fit. Once 
again fairly good agreement is obtained throughout the angle-of -attack 
range . 

The same comparison is made in figure 11 as in figure 10 except that 
the 1-3-5 approximation of reference 3 and an accompanying least squares 
fit have been used. Although relatively more scatter was found in the 
experimental data when reduced by the method of reference 3, fairly good 
agreement between the two methods is noted in this figure. Percentagewise, 
this comparison indicates about 10 percent difference in the moment coeffi- 
cient at the highest angle of attack. This is about the same percentage, 
and in the same direction, that Murphy' s solution for the period of the 
oscillation differs from the exact solution for the corresponding case 
of planar motion governed by a 1-3-5 moment. 

It is perhaps desirable at this time to summarize briefly the 
advantages and disadvantages of the three methods that have been 
considered. 

1. Rasmussen's method (ref. 2) has the advantages of not being 
restricted to planar (or reasonably planar) motion and of being an exact 
solution. It has the disadvantages of being restricted to the 1-3 case 
only and of being relatively difficult to apply except in its inexact 
form. 


2. Murphy's method (ref. 3) has the advantages of not being 
restricted to planar (or reasonably planar) motion, of having the widest 
range of permissible nonlinear moment approximations, and of being rela- 
tively easy to apply. It has the disadvantage of being an approximate 
method from a mathematical viewpoint, with its range of applicability not 
clearly defined. 

3. The method of the present report has the advantages of being an 
exact solution, of having a fairly wide range of permissible nonlinear 
moment approximations, and of being relatively easy to apply. It has the 
disadvantage of being restricted to reasonably planar motion. 


CONCLUSIONS 


A method has been presented for obtaining the aerodynamic stability 
characteristics of bodies of revolution that are governed by nonlinear 
restoring moments from data obtained from free-flight tests. Four dif- 
ferent approximations have been presented and solved in closed form. The 
pertinent conclusions that result from this analysis are as follows: 
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1. Although the derivations are based on the existence of planar 
motion, with zero roll rate and damping, the solutions have been found 
to be applicable when the roll rate and damping are small and when the 
motion in the a, 0 plane can be classed as reasonably planar. Under 
these conditions, the method yields realistic values of the moment 
coefficient as a function of angle of attack. 

2. The various approximations can be programmed for rapid solution 
on a digital computer with only a modest amount of input information 
necessary. 

3. The resulting solutions are, to a certain extend, self-checking, 
since only two (or three) data points are used to obtain a given approxi- 
mation and additional data points can be used to verify the adequacy of 
the approximation. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Jan. 12, 1961 
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Figure 1.- Effective angle-of-attack parameter for a single term nonlinear restoring moment. 
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Figure 2.- Effective angle-of-attack parameter for a linear plus cubic restoring moment. 
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Figure 3 • - Cases existing under a linear plus cubic restoring moment assumption. 
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Figure 4.- Cases existing under a 1-3-5 and a 1-2-3 restoring moment assumption. 
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Figure 5«- Moment-curve slope from linear theory as a function of the 
maximum resultant angle of attack. 
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Figure 6.- Comparison of 1-3-5 &ncL 1-2-3 fits of experimental data points. 
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Figure 8.- Comparison of free-f light data with wind-tunnel data. 
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Figure Comparison of the method of the present report with Rasmussen's 

method . 


46 



Figure 10.- Comparison of the method of the present report with Murphy's 

1-3 method. 
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Figure 11.- Comparison of the method of the present report with Murphy's 

1-3-5 method. 
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